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We investigate the growth of the total number of particles in a symmetric exclusion process driven 
by a localized source. The average total number of particles entering an initially empty system grows 
with time as yt in one dimension, i/lni in two dimensions, and linearly in higher dimensions. In 
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are independent of the intensity of the source. We also discuss fluctuations of the total number of 
particles and determine the asymptotic growth of the variance in one dimension. 
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I. INTRODUCTION 



average (N) grows as 



A system of particles undergoing symmetric random 
walks on a lattice under the constraint that multiple oc- 
cupancy is forbidden, a so-called symmetric exclusion 
process (SEP), has been extensively investigated, espe- 
cially in one spatial dimension (see, e.g., books and re- 
views [THZ])- A surprising feature of the SEP, at least at 
first sight, is that the average density evolves as if there 
were no exclusion. (More subtle characteristics are af- 
fected by exclusion, e.g., the mean-square displacement 
of the tagged particle exhibits a remarkable sub-diffusive 
growth, (x 2 ) ~ Vt) m one dimension [SHU]; in higher 
dimensions, the diffusion behavior of the tagged particle 
is restored and only the amplitude of the mean-square 
displacement acquires the dependence on the concentra- 
tion.) This may lead to the impression that the exclu- 
sion property is irrelevant as long as the average charac- 
teristics are concerned. Further, the peculiarities of the 
SEP arise in one dimension; when d > 1, the SEP and a 
system of non-interacting random walks exhibit similar 
qualitative behaviors. 

The above conclusions hold in the absence of sources. 
If particles enter the system through external reservoirs, 
the exclusion property may become crucial. Here we ana- 
lyze the SEP with a localized source and find that exclu- 
sion plays an important role, particularly in low dimen- 
sions d < 2. More precisely, we show that the critical 
dimension is d c — 2, viz. similar behaviors arise when 
d > 3. (When d — 2, the difference with the higher- 
dimensional behavior is logarithmic, i.e., rather small.) 

We study the SEP on the d— dimensional hyper-cubic 
lattice Z d . We set to 1/d the hopping rate to each of 
the 2d neighboring sites, so that the total hopping rate is 
equal to 2 for all d. (The hopping event is allowed only 
when the destination site is empty.) We denote by F the 
flux of particles to the origin, that is, the rate at which 
new particles would be arriving if the origin were always 
empty. We want to understand how the total number of 
particles N(t) grows. (The source is turned on at t = 0; 
we assume that the system is initially empty.) 

The quantity N(i) is a random variable. In the large 
time limit, N(t) is concentrated near its average, that 
is, fluctuations are relatively small. We show that the 
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in the long time limit. Intriguingly, the leading asymp- 
totic behaviors in one and two dimensions are indepen- 
dent of the flux F as long as it is positive. Also, in one 
and two dimensions the total number of particles enter- 
ing the system is a negligible fraction of the number of 
particles, Ft on average, which would have entered in 
the no-exclusion case. In three and higher dimensions, 
the average total number of particles entering the system 
grows linearly with time. The renormalized flux <f>d(F) 
admits a simple expression via the 'bare' flux F and a 



Watson integral Wd [see ( 12 )] 
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Thus the renormalized flux grows linearly with bare F 
when F -C 1, and saturates to W^ 1 in the F — > oo limit; 
needless to say, & d (F) < F. 

We also discuss fluctuations of the total number of in- 
jected particles. Specifically, we argue that the variance 
(N 2 ) c = (N 2 ) — (N) 2 exhibits an asymptotic growth 
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and we determine the amplitude V% in one dimension 
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The computation leading to this result is performed in 
the F = 00 limit, but we argue that the amplitude is 
independent on the flux F. 

We study the SEP with a localized source, but the 
results are more widely applicable. For instance, some 
stochastic processes are mathematically similar, or even 
isomorphic, to the SEP with a localized source with 
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infinitely strong flux. Examples include monomer- 
monomer catalysis and the voter model [T3lJ15| . Both 
models are particularly natural in two dimensions and 
they have been analyzed for d < 2; some of the results of 
Refs. [TSTfTB"] are equivalent to the predictions in one 
and two dimensions. The spreading of very thin wetting 
films has been described by a SEP-like model (see [IB] for 
a review) . The natural dimensionality of the substrate is 
d = 2 and the average injected mass has been computed 
[17] for d < 2 , and has been shown to compare favorably 
with experimental observations. More complicated SEP- 
like models with a localized source have also been stud- 
ied. One recent example [TS] involves synthetic molecu- 
lar motors, so-called molecular spiders [19]; in contrast 
to particles in the SEP, the underlying model of molecu- 
lar spiders |20j is non-Markovian as the motion of spiders 
affect the substrate. 

The rest of this paper is organized as follows. A heuris- 
tic explanation of asymptotic behaviors is given in 
Sec. [TT] In Sec. |III| we present an exact solution in the 
limiting case of the infinitely strong source, from which 
we determine exact asymptotic behaviors. The general 
case of finite flux is treated in Sec. IIVI where we derive 
the renormalized flux, Eq. ([2]), when d > 2, and we com- 
pute sub-leading corrections in one dimension. In Sec. [V] 
we discuss fluctuations. In particular, we analytically 
determine the variance in the one-dimensional setting. 
Section |VI| contains a summary. 



II. HEURISTIC DERIVATION 

Here we give a heuristic explanation of the growth laws 
0. We consider the simplest situation when the source 
is infinitely strong, F = oo. In this case, the density at 
the origin is po = 1 for all t > 0. The initial density in 
other sites is zero. In one dimension, the density decays 
from 1 near the origin to zero on distances exceeding the 
diffusive scale y/i. The density profile is roughly linear 
in the diffusive layer: 



l-\x\/y/i \x\ < Vt 
\x\ > Vi 



(5) 



Therefore, (N) ~ in agreement with ([I]). 

In two dimensions, the density is also time-dependent. 
Using the well-known asymptotic (see, e.g., |21| ) 



P ' 



lnr 



(6) 



we obtain 
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In three dimensions, the density becomes stationary. 
Thus one must solve the Laplace equation V 2 p = 0. An 
appropriate solution is p(r) ~ r _1 (Coulomb's potential), 



or more generally p(r) ~ r -( d - 2 ) when d > 2. There are 
no particles on distances far exceeding \ft. The average 
number of particles is estimated by integrating over the 
ball of radius of the order of \fi. This gives 
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where in the last step we used p(r) 
growth law ^ applies when d > 2 
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III. INFINITELY STRONG SOURCE 

The case of an infinitely strong source is exactly solv- 
able as far as the average quantities are concerned. The 
average density p x (t) satisfies the diffusion equation 
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Here x = (x\, . . . , %d) € Z d is the lattice site and V 2 de- 
notes the discrete Laplace operator; e.g., in one dimen- 
sion V 2 pa; = Px+i — %Px + px—i- The entire hopping rate 
is set to 2, so the individual hopping rates are equal to 
1/d as any site of the hyper-cubic lattice has 2d nearest 
neighbors. The initial condition reads 



p x (t = 0) = 5^o 
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An infinitely strong flux is represented by the boundary 
condition 
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assuring that there is always a particle at the origin. 
Equations ([8]) — (p~0|) admit an exact solution, from which 
we will extract leading asymptotic behaviors of the aver- 
age total number of particles: 



d=l, 
2irt/\nt d = 2, 
t/W d d>2. 

Here Wd are the so-called Watson's integrals [22] 
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where q = (qi, . . . , qj) and 
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Watson's integrals often appear in problems involving 
lattice Laplacians. In three dimensions, the Watson in- 



tegral (12 1 has been expressed [53] via Euler's gamma 
functions: 
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The solution to (|8])-( 10 1 can be established using var- 



ious approaches. Here we outline a derivation following 
an approach of Refs. [T2HT5] which holds in arbitrary spa- 
tial dimension and gives exact results for Laplace trans- 
forms of the basic quantities. Asymptotic behaviors of 
the Laplace transforms imply asymptotic large time be- 
haviors. We then use another method and obtain an 
explicit solution in one dimension in a more direct way, 
namely, without using the Laplace transform. 



A. Solution via Laplace Transform 

A solution to Eqs. ([8])-([9]), i.e., the solution of the prob- 
lem without the boundary condition (10), is merely the 
lattice Green function I x (2t/d)e~ 2t . Here 



IP- 
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(14) 



is the shorthand notation for the product of the mod- 
ified Bessel functions. To maintain the validity of the 
boundary condition po(t) = 1 throughout the evolution, 
we employ a simple trick: We add particles to the ori- 
gin at a certain rate & d (t) which we choose to fulfill the 
boundary condition. Since the governing equation is lin- 
ear, the general solution to ([8])-([9]) with a source is a 
linear combination 
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The integral term is the contribution due to the source 
<& d (r)dT which is added at the origin during the time 
interval (t, t + dr); the original source S(t) yields the 
first term on the right-hand side of (15). At the origin 
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This integral equation determines the strength § d (t) of 
the source. To extract an explicit expression we notice 



the convolution structure of the integral in (16). Hence 



we apply the Laplace transform. This yields a neat rela- 
tion 
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between the Laplace transform of the strength of the 
source and the Laplace transform of the power of the 
Bessel function Iq: 
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B d (s) = I dte- st [l (^)] d e- 2t (18b) 



Using the integral representation of the modified Bessel 
function, Io(t) = f ™ ^ exp(tcosg), and performing the 
integration over t, we re-write ( 18b ) as an integral 
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with Q(q) given by (13). Expanding the right-hand side 



of (j 1 9[) we determine the small s behavior 
-l 

Bd(s) 
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where B d (0) = W d for d > 2. Combining Eqs. (17) and 



( 20 1 we find the leading s — > behavior of the Laplace 



transform of the strength of the source 

Ms) - 
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(Was)- 1 d>2, 



Converting (21) we obtain the large time asymptotic of 



the strength of the source 
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By inserting ([22]) into the relation (N(t)} = / dr$ d (r) 
we arrive at the announced asymptotic (11). 



Using Eqs. (15) and (22) wc can additionally extract 



the asymptotic behavior of the density. In one dimension, 
it is actually simpler to employ a continuum approach 
from the very beginning. Thus we must solve the initial- 
boundary value problem 

d t p = d xx p, p(x = Q,t) = l, p(x,t = 0) = 

for the average density p(x,t). The solution is 



where erfc(u) 



p{x, t) — erfc 
f°° dv e _1 
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1 — erf(it) is an er- 
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ror function. The asymptotically exact result ([23]) is the 
corrected version of (5|. 

In deducing the asymptotic behavior of the density in 
two dimensions, we first recall that 
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We now simplify the dominant integral term in Eq. ( 15 ) 
by using (l22j) and plh to yield 



p x (t) = (In*)" 1 E 1 (£ t ) 
x|.and E 1 (z) = ^ d f 
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Here r 2 = x\ + x\ and E\{z) = — e _zn is an expo- 
nential integral. Using the asymptotic behavior of the 
exponential integral, E\{z) = — In2 + G(l) as z — > [24] . 
we see that Eq. ([5| provides the leading asymptotic not 
too far from the source, namely when r -C yt. 

In three and higher dimensions it suffices to consider 
the final (stationary) density. Therefore we put t = oo 
into Eq. (Il5| and find 
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In deriving Eq. (27) from Eq. (|15j we have also used the 



integral representation 
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dq 
2^ 
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of the Bessel function. For d > 2, the steady state solu- 
tion (27 1 is clearly non-trivial. For instance, at the sites 



closest to the origin the density is p' = 1 — (2W<j) , in 
three dimensions, p' w 0.34053733. 



B. Explicit Solution in One Dimension 

Due to the p-j = pj symmetry, it suffices to consider 
the densities pj to the right of the origin, j > 0. Thus 
we must solve 



Pj = Pj+i - 2 Pj + Pi- 



ty 



(28) 
and the 



subject to the initial condition pj(t 
boundary condition ( |10[ ). We can instead consider (28) 
on the entire lattice and impose the initial condition 



Pj(t = 0) = { 



j<0 
3=0 
j>0 



(29) 



Then the boundary condition po — 1 will manifestly hold. 
(The resulting solution is applicable, of course, only to 
j > 0; for j < 0, proper densities are recovered from 
relation p_j = pj.) The solution of the discrete in space 
diffusion equation (28) subject to the initial condition 



(29) is straightforward: 
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The average total number of particles to the right of the 
origin is J2j>oPji an d therefore 



(N) = p + 2j2pi 

j>0 



(31) 



Plugging (30) into (31) and massaging the sums we get 
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The sum on the right-hand side of Eq. ( 32 ) can be ex- 



pressed through the Bessel functions Io and 1\ to give 

(N) = e- 2t [I„(2t) + 4il (2i) + 4^(2*)] (33) 

The leading large time behavior of (N(t)) agrees with 
( 11 ), and as a byproduct of having a compact exact result 



(33) we can also extract sub-leading corrections: 
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IV. FINITE FLUX 

Here we investigate the general case when the flux is 
finite, F < oo. In one dimension, for instance, we need 



to solve Eq. ( 28 1 away from the origin, while at the origin 



the average local density obeys 



p = 2(p 1 -p )+F(l-po). 



(34) 



Due to the symmetry, pj = p~j, and it suffices again to 
consider only j > 0. 

In one and two dimensions the leading asymptotic be- 
havior is independent of F for any F > 0. For instance, 
in one dimension p — > and {p\ — po) — > as t — > oo, 
and therefore Eq. (34) tells us that po 1. Thus the 
boundary condition (10) is asymptotically correct, and 



hence the leading asymptotic behaviors of all densities 
and of the total average number of particles are the same 
as in the case of the infinitely strong source. The two- 
dimensional case is treated similarly. In three and higher 
dimensions, the average density at the origin exceeds the 
average densities at neighboring sites even as t — > oo and 
hence po < 1. 



A. Asymptotic Behavior in One Dimension 

To investigate the emergence of the asymptotic behav- 
ior in a quantitative manner we again employ the Laplace 
transform. Equations (28) and (34) become 



spj = p j+1 - 2pj + pj_x , j > 
sp = 2(pi - po) + Fis^ 1 - po) 



(35) 



where 



Pj(s) = 
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Equations ( 35 ) admit an exponential solution 
p^s) = A(s)[a(s)Y 



(36) 
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Plugging this ansatz into ( 35 1 we get 

s + 2 - Vs 2 + 4s 
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The total average number of particle is 

oo oo 

(N(t))= Pj(*) = "/*>(*) + 2 

3=—oo j=0 

and therefore its Laplace transform is 



(37) 



t (N(t))=A(s) 
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l + a(s) 
1 - a(s) 
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Using (371 we expand the right-hand side of (381 in the 
s — > Q limit to yield 
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Using this expansion in conjunction with (38) we deduce 



the large time expansion of the average total number of 
particles 



(N(t)) 
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This asymptotic expansion confirms our assertion that 
the leading behavior is universal (that is, independent 
on F). 

We can also find the large time behavior of the densi- 
ties to substantiate the claims po — ¥ 1 and p\ — > 1 which 
have been previously made. Indeed, using (|36[) and (37) 



we establish the small s expansions of the Laplace trans- 
forms 
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from which we deduce the t — > oo behaviors 
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B. Renormalized Flux in Three and Higher 
Dimensions 



In three and higher dimensions, the flux F affects lead- 
ing asymptotic behaviors. The renormalized flux is time- 
independent in the long time limit, and therefore (15) 
gives 



Px = <S>d I dTl x (^)e~ 
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(41) 



in the long time limit. Recalling (14) we see that the 
density at the origin is 

Po = / dr [I (^)] d e- 2r = 3> d W d (42) 
Jo 

On the other hand, the renormalized flux is related with 
the density at the origin via 



$ rf = F(l-p ) 



(43) 



Combining ( |42| ) and ( 43 1 we arrive at the announced ex- 
pression ([2]) for the renormalized flux. 



FLUCTUATIONS 



The average numbers of injected particles [see Eq. ( 11 )] 



are lattice- independent when d < 2. This occurs because 
in one and two dimensions the spatial scale where the 
density varies is growing with time, so that the lattice 
structure is asymptotically irrelevant. Hence the average 
quantities, e.g. the average densities [Eqs. (23) and (26)], 
can be established in the realm of continuum approaches. 
To probe fluctuations in low dimensions, d < 2, one can 
also utilize continuum approaches. 

For d > 3, the results are lattice-dependent. The av- 
erage density is stationary when d > 3 and this feature 
simplifies the problem. Namely, the total number of in- 
jected particles should be a Poisson distributed random 
quantity. In the long time limit, it becomes a Gaus- 
sian distribution with equal (in the leading order) aver- 
age and variance (N 2 ) c = (N 2 ) — (N) 2 . In one and two 
dimensions, one still expects the average (N) and the 
variance (iV 2 ) c to exhibit the same dependence on time, 
but (N) ^ (iV 2 ) c as there is no reason for the equality. 

These arguments lead to the announced results p]). 
The amplitudes V\ and V2 characterizing the variance in 
one and two dimensions should not depend on the flux 
F, The reason is the same as in the case of the average 
(N) , although for the latter quantity the evidence is much 
stronger, and in one dimension we even have an exact 
solution (for the Laplace transform) and the asymptotic 
expansion ( 39 ) which explicitly show that F affects only 



sub-leading corrections. 

In this section we derive the amplitude in one dimen- 
sion, Eq. Q. In the computation we set F = 00; as 
we stated above, the result Q is apparently universal, 
namely it holds for any F > 0. 

The variance can be computed by various techniques. 
In one dimension, one can use exact methods. To deter- 
mine the most interesting long time behavior, it is easier 
to employ continuum approaches. These methods are, 
in principle, applicable both in one dimension and two 
dimensions, but the former case is much more tractable, 
so we consider only a one-dimensional setting. We shall 
compute the variance using fluctuating hydrodynamics 
(see [TJ [5S]). According to this approach one should 
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solve a Langevin equation [TJ [5] 

d t q = d xx q + d x y / a(q)^(x,t) 



(44) 



for the fluctuating particle density q(x, t). Here £(x, t) is 
a Gaussian white noise with standard correlations 

{Z(x,t))=0, ^(x,t)^x',t'))=S(x-x')6(t-lf) (45) 

The quantity u{q) characterizing fluctuations of the cur- 
rent is known [II [5] to be a(q) = 2q(l - q) for the SEP. 

When F = oo, there is no interaction (for all t > 0) 
between regions to the left and to the right of the source. 
Consider the x > region. Linearizing the Langevin 
equation ( 44 1 around the hydrodynamic solution ( 23 1 , 



i.e., writing q = p + qi and assuming that q\ <C p, we find 
that the perturbation obeys a diffusion equation with a 
stochastic source 



d t qi - d xx qi = d x y/<r(p) £(x, t) 



Setting qi = d x ip, we rewrite (46) as 



d t ip - d xx ip = y/ a(p) £(x, t) 



(46) 



(47) 



We must solve this equation subject to the initial con- 
dition ip(x,t = 0) = and the boundary condition 
d x ij){x = 0,t) = which follows from the requirement 
q(x = 0, t) = p(x = 0, t) = 1. The solution reads 



Jo Jo v 4 »C - f) 
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4(t-t') _|_ g 4(t-t') 



(48) 



The total number of particles N+ in the x > region is 

poo 

N + = (N+)+ / dx qi (x,t) (49) 
Jo 

The variance is therefore 



(N 



Do 



dx / dyqi(x,t)q 1 (y,t) 
o Jo 



dx / dyd x ^(x,t)d v ^(y,t) 



= <V< 2 (o,t)> 



(50) 



Plugging ( 48 ) into ( 50 1 and using ( 45 ) we obtain 
rt dt' 



W)c = 



1 
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dya[p(y,t')]e 2 <*-'') (51) 



The total number of particles iV_ in the x < region 
is a random quantity which is independent on 2V+ and 
identically distributed. Therefore for the total number of 
particles N = N- + N + the variance is (N 2 ) c - 2{N^) C . 
Using this together with (51 1 we obtain 



(N 2 ) c = 



*dT 



dx cr[p(x, t — t)] e 



(52) 



Making the transformation t — tT and x — y/tX, and 
using a(p) = 2p(l — p) together with Eq. (23 1, we recover 
the time dependence (iV 2 ) c = V\\ft with amplitude 



V 1 = ±{U 1 -U i ) 



(53) 



where 



dT 



dX 



erf 



X 



Computing the first two integrals, 

Ui = v 7 ^ - x/^ , U 2 = (3-2V2)v / 2^ : 
we simplify (53 1 to the announced result Q. 



VI. SUMMARY 

We studied the growth of the total number of particles 
in a symmetric exclusion process driven by a localized 
source. Specifically, we assumed that new particles are 
injected into a single lattice site, the origin, whenever 
the origin is empty. We showed that the average total 
number of particles entering an initially empty system 
exhibits a simple asymptotic growth 0. In one and two 
dimensions, the leading asymptotic behaviors for the av- 
erage total number of particles turn out to be universal 
(independent of the flux) and in both cases the asymp- 
totic growth is slower than linear in time. In three and 
higher dimensions, the average total number of particles 
entering the system grows linearly with time, namely as 
<&d(F) t. We derived a simple equation ^ expressing the 
renormalized flux 3>d(F) through the bare flux F. 

In one and two dimensions, the results are insensitive 
not merely to the flux F, but also to the detailed struc- 
ture of the source (the injection may occur through a few 
sites; the number of such sites and the details of their lo- 
cation are irrelevant) and to the lattice structure (in two 
dimensions, there is no need to require that the lattice is a 
square grid). In three and higher dimensions, the results 
are sensitive to the aforementioned detailed properties. 
It seems that if the lattice is arbitrary, but new particles 
are still injected into a single lattice site, the expression 
([2]) for the renormalized flux remains valid, although one 
has to use a proper expression for the Watson integral 
depending on the underlying lattice; e.g., for the body 
centered cubic lattice (bcc) the Watson integral is 



d 3 q 



w 2(2-7r) 3 1 — cosgi cos (72 cos ^3 
\T(-)] 4 

1 u;j = 0.696601966... 



(2tt) 2 



We also discussed fluctuations of the total number 
of particles, especially the variance for which we gave 
asymptotic growth laws ([3]). The arguments leading to 
([3]) are solid only in one dimension where we used fluc- 
tuating hydrodynamics and analytically established the 
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asymptotic growth of the variance. In two dimensions, 
the result cited in Eq. ^ is just a guess with an unknown 
multiplicative factor V" 2 i even the functional form of the 
variance is conjectural. It would be interesting to com- 
pute the variance in two dimensions by employing fluc- 
tuating hydrodynamics or another continuum approach 
[25j . In three and higher dimensions, the prediction of 
Eq. g for the variance is sharp yet unproven; one would 
like to justify that prediction, or disprove it. 

The one-dimensional case is particularly tractable, and 
one may be able to determine higher cumulants. To 
guess the outcome, we notice a similarity of our prob- 
lem and the problem of the evolution of the SEP start- 
ing with a step-function initial condition, particularly 
p(x, t = 0) = 1 for x < and p(x, t = 0) = for x > 0. 
In this latter problem all cumulants of the total current 
grow in a diffusive manner 26, 27J. The same should be 
valid for the cumulants of N(t) in our problem. Since 
the numbers N + and JV_ of particles to the right and left 
of the origin are independent (in the simplest case of in- 
finite flux) and identically distributed random variables, 
it suffices to consider N + (t). Similarly to the SEP with 



a step-function initial condition |26) we anticipate that 



P(N+,t) = Prob 



y/i 



-y/tG{n) 



(54) 



The derivation of the large deviation function for the 
SEP with the step-function initial condition |26) is com- 
plicated and the results do not seem to admit a straight- 
forward extension to the present case. For instance, the 
probability to have zero total current is oc e _v/ * |26| . 
which means that the large deviation function remains 
finite in the small current limit: G s t cp (0) < oo. In our 
case, N + (t) = 0, or more generally that N + (t) is small, 
with probability oc e _t ; this implies that the large de- 
viation function diverges, G(n) ~ n" 1 as n — > 0. Both 
settings are still very close, so the methods of Ref. [26] 
could be applicable to the present case, at least in the 
situation when the flux is infinite. 

I thank Kiron Mallick, Baruch Meerson, Gleb Oshanin, 
and Darko Stefanovic for useful correspondence. 
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